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Multi-oscillatory beliavior of mitochondrial inner membrane potential in self- 

organized cardiac mitochondrial networks can be triggered by metabolic or oxidative 
stress. Spatio-temporal analyses of cardiac mitochondrial networks have shown that 
mitochondria are heterogeneously organized in synchronously oscillating clusters in which 
the mean cluster frequency and size are inversely correlated, thus suggesting a modulation 
of cluster frequency through local inter-mitochondrial coupling. In this study, we propose 
a method to examine the mitochondrial network's topology through quantification of its 
dynamic local clustering coefficients. Individual mitochondrial A*rn oscillation signals were 
identified for each cardiac myocyte and cross-correlated with all network mitochondria 
using previously described methods (Kurz et al., 2010a). Time-varying inter-mitochondrial 
connectivity, defined for mitochondria in the whole network whose signals are at least 
90% correlated at any given time point, allowed considering functional local clustering 
coefficients. It is shown that mitochondrial clustering in isolated cardiac myocytes changes 
dynamically and is significantly higher than for random mitochondrial networks that are 
constructed using the Erdos-Renyi model based on the same sets of vertices. The 
network's time-averaged clustering coefficient for cardiac myocytes was found to be 
0.500 ± 0.051 (A/ = 9) vs. 0.061 ± 0.020 for random networks, respectively. Our results 
demonstrate that cardiac mitochondria constitute a network with dynamically connected 
constituents whose topological organization is prone to clustering. Cluster partitioning in 
networks of coupled oscillators has been observed in scale-free and chaotic systems 
and is therefore in good agreement with previous models of cardiac mitochondrial 
networks. 
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INTRODUCTION 

The mitochondrial network in cardiac myocytes consists of highly 
organized and densely packed mitochondria whose inner mem- 
brane potential A>I'„, can be triggered to enter several cycles of 
de-and re-polarizations by numerous stressors such as oxida- 
tive or metabolic stress (see Aon et al., 2008a for a review). 
These oscillations can be strictly localized in the form of transient 
single mitochondrial depolarizations (Nivala et al., 2011), indi- 
vidual or clustered mitochondrial A^m oscillations (Romashko 
et al, 1998; Kurz et al., 2010a) with clusters that can span the 
whole myocyte (Aon et al., 2004). Recruiting neighboring net- 
work mitochondria into an initial synchronized nucleus of a few 
mitochondrial oscillators has been described to be a fundamen- 
tal process for global network synchronization (Strogatz, 2000, 
2003; Aon et al, 2008b). During this process, a mitochondrial 
cluster can reach a critical size (sometimes referred to as "mito- 
chondrial criticality" Aon et al., 2004, 2006) where mitochondria 
spontaneously self-synchronize, as in a phase transition. So far, 
investigations strongly support the fact that ROS-induced ROS 
release is a key player in such inter-mitochondrial communica- 
tion or coupling (Zorov et al., 2006; Zhou et al., 2010; Nivala et al., 
2011). 



Recently, wavelet-based analysis tools have been developed 
to examine the mitochondrial network's spatio-temporal behav- 
ior under pathophysiological conditions (Kurz et al., 2010a,b); 
dynamic frequencies could be allocated to individual mitochon- 
dria and clusters of mitochondria with similar frequencies were 
identified that allowed for a quantitative characterization of 
the cluster's network properties. However, the network's topol- 
ogy (or connectivity properties) as opposed to its architec- 
tural organization (see Aon and Cortassa, 2012 for a review) 
has not yet been investigated quantitatively in terms of its 
clustering properties. Mitochondria in cardiac myocytes serve 
as the main energy supplier and modulator of the myocyte's 
mechanical and electrical processes, but are also modulated 
by the latter; therefore, the mitochondrial network's topolog- 
ical heterarchy becomes increasingly complex and non-linear 
(Yates, 1993). The functionality of an individual mitochondrial 
network node, though, can in part be characterized through 
its connectedness with other network nodes, (cf. Passingham 
et al., 2002), a result of the interplay of the entire complex 
mitochondrial network as an integrated system. The clustering 
coefficient can be used as a measure of the network's robust- 
ness toward the functional deletion of single mitochondria or 
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the network's efficiency to communicate (metabolic or other) 
information. 

The present work investigates the presence of functional 
(dynamical) connectedness in the form of clustering of mito- 
chondrial networks in isolated cardiac myocytes in comparison 
with clustering in random networks based on the Erdos-Renyi 
model. Functional network clustering is subsequently related to 
the network's spatio-temporal properties of the major cluster of 
mitochondria with similar frequencies. 

MATERIALS AND METHODS 
EXPERIMENTAL METHODS 

All experiments were carried out on freshly isolated adult 
guinea pig ventricular myocytes at 37° C following protocols that 
were previously described (O'Rourke and Marban, 1994) with 
approval from the Johns Hopkins University Animal Care and 
Use Committee and in accordance with guidelines established in 
the Guide for the Care and Use of Laboratory Animals, published 
by the National Institutes of Health (NIH Publication No. 85-23, 
Revised 1996). 

In brief, cardiac myocytes were perfused with Tyrode solu- 
tion (pH 7.5) containing 1 mM Ca^+ in the presence of 10 mM 
glucose and oscillations were triggered with a localized (5 x 
5 |xm) laser flash, as previously described (Aon et al., 2003). 
Mitochondrial inner membrane potential was monitored 

with the cationic potentiometric fluorescent dye tetramethylrho- 
damine methyl ester (TMRE) and images were recorded with a 
two-photon laser-scanning microscope (MRC-1024MP, Bio-Rad) 
with excitation at 740 nm (Tsunami Ti:Sa laser, Spectra-Physics) 
and red emission of TMRE was collected at 605 nm using a band 
pass filter 578-630 nm (Aon et al, 2003). 

SELECTION AND PROCESSING OF INDIVIDUAL MITOCHONDRIAL TMRE 
SIGNALS 

As detailed before (Kurz et al, 2010a,b), TMRE signals of individ- 
ual mitochondria were extracted from planar images of isolated 
cardiac myocytes recorded at a rate of dt by manually applying a 
grid template on a pixel-by-pixel basis to the averaged image of n 
subsequent images in time starting at the onset of TMRE oscilla- 
tions. Parameter n was chosen such that n ■ dt was smaller or equal 
to the smallest period of all TMRE oscillations. Only myocytes 
with no shifts in the z-direction were included and shifts in the 
X and y-directions were corrected by moving the template grid 
accordingly. 

The subsequent spatio-temporal signal analysis was conducted 
using the wavelet transform, correlation and coherence analysis 
as previously described in detail (Kurz et al, 2010a) and outlined 
below. 

WAVELET ANALYSIS 

With no prior knowledge of whether mitochondrial oscillations 
are stationary, wavelet analysis has been used to probe the dynam- 
ically changing frequencies of cardiac mitochondria (Kurz et al., 
2010a). In this study, the same wavelet analysis was utilized: for 
each mitochondrion's TMRE signal, Morlet wavelets were taken 
where spacing between scales was set to dj = 0.1 and the small- 
est scale of the wavelets was chosen as sq = 4df, signifying the 



smallest possible period for the detection of one oscillation. The 
total number of scales was determined as ji = \0g2iN / sq) / dj + 1 
with N being the total number of the recorded images per cell, 
thus resulting in scales that range from sq to sqI^^^^'''^' and each 
scale having dj suboctaves. To determine cutoff frequencies, the 
longest period T of a synchronized oscillation of one cell was 
identified to determine the minimum cutoff frequency as Vmin = 
1 / 1 . 1 r and, similarly, the maximum cutoff frequency at Vmax = 
I/sq. Power lineplots between Vmin and Vmax were interpolated 
for every scale with segments of 0.1 mHz and the maximum 
power for the interpolated plots was determined to eventually 
obtain maximal scale frequencies at each time point for each 
mitochondrion. 

In Figure lA, three random mitochondria labeled 1, 2, and 3 
from an isolated cardiac myocyte are chosen to illustrate individ- 
ual mitochondrial TMRE signal behavior. Evidently, only mito- 
chondria 1 and 3 show marked oscillatory behavior over time 
whereas mitochondrion 2 is non-oscillating (see Figures 1B,C). 
Figure IC shows the absolute squared wavelet transform of mito- 
chondrion 3 over frequency and time. It can be seen that the 
main frequency component, depicted by the dark red color, varies 
between 15 and 20 mHz during the recording. This frequency 
component corresponds to the time interval between oscillation 
peaks and troughs observed in the upper panel of Figure IC, that 
is approximately between 50 and 65 s. 

SELECTION OF MITOCHONDRIA BELONGING TO A MAJOR CLUSTER 

By obtaining frequency histograms of all mitochondria for each 
time point t, the maximum peak of each histogram was identi- 
fied and other peaks were considered significant if their histogram 
amplitude was above 10% of the respective maximum peak (see 
also Figure 2A in Kurz et al., 2010a). If the mean TMRE signal 
of all mitochondria corresponding to the maximum peak had a 
cross correlation of 95% or higher with the mean TMRE signal 
of mitochondria in an adjacent peak, the latter was incorpo- 
rated in the maximum peak and the procedure repeated with the 
next adjacent peak until the correlation dropped below 95%. The 
respective signal cross correlation was determined over a running 
window T„ = I. IT with center at time-point f . In the resulting 
maximum peak at each f, all mitochondria in that cluster were 
identified and their mean TMRE signal in T-^, was cross-correlated 
with that of each mitochondrion in the myocyte that did not 
belong to the cluster. Again, if the correlation coefficient was 95% 
or higher, the individual mitochondria were incorporated into 
the cluster, thus yielding the major cluster of mitochondria that 
are highly correlated at time t. The normalized dynamic area of 
the major cluster was taken as the quotient of the total dynamic 
cluster pixel count and the total myocyte pixel-count. 

COHERENCE ANALYSIS 

Temporal properties of major cluster mitochondria were exam- 
ined by considering the average coherence of each cluster mito- 
chondrion's TMRE signal with all of its nearest neighbors. 
Coherence values range between one and zero, indicating whether 
two signals oscillate synchronously at each frequency or not, 
respectively. Using a running window T„, a fixed Fast Fourier 
Transform length of 2^^ dt, a frequency range for each myocyte of 
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FIGURE 1 I Single mitochondrial signals and wavelet analysis. (A) 

Mitocliondria in cardiac myocytes are densely packed and their individual 
TMRE signal can be extracted as previously described (Kurz et al., 2010a, b). 
(B) Two mitochondrial signals with different oscillatory patterns from 
different locations within the myocyte are shown. (C) Absolute squared 
wavelet transform over frequency and time of an oscillating mitochondrion. 
The major frequency component varies between 15 and 20mHz. 



0-100 mHz and its division into {2^^ /2) + 1 segments (each seg- 
ment tlierefore corresponding to ^O.l mHz), the coherence was 
determined between each mitochondrion and its nearest neigh- 
bors for each time point and the mean time-resolved coherence 
over all nearest neighbors was calculated. Consequently, the aver- 
age coherence for each major mitochondrial cluster at each time 
point was obtained. To compare myocytes with unequal duration 
of recordings, the duration of oscillation of each cardiac cell was 
set to 1. 



MITOCHONDRIAL CLUSTER AMPLITUDE 

For the mean TMRE signal of the major cluster, aU peaks and 
troughs were manually identified and the respective peak-trough 
amplitudes and peak-trough time differences were determined 
(see also Figure 4 in Kurz et al., 2010a). Normalization of peak- 
trough amplitudes was achieved by dividing the amplitudes 
through the respective maximum amplitude. 

LOCAL CLUSTERING COEFFICIENT 

In graph theory, clustering describes the functional topology 
of a network by quantifying the degree to which a set of net- 
work vertices (or constituents) resembles the graph-theoretical 
concept of a clique (i.e., a set of vertices in which every two 
vertices are connected to each other). Considering a single net- 
work vertex m, the fraction of those network vertices that are 
(undirectedly) connected to vertex m, with respect to the maxi- 
mum number of topological neighbors of vertex m, are described 
by the clustering coefficient Evidently, if vertex m has nip] 
topological neighbors, then the maximum number of possible 
undirected links among these neighbors is misiitnN — l)/2. With 
L„, being the number of undirected links between the neighbors 
of mitochondrion m, local clustering Cm is hence defined as Cm = 
2Lm/ (misiimis! — 1)). According to Watts and Strogatz (1998), the 
mean clustering coefficient of the whole complex network C is 
determined by the arithmetic mean of all local clustering coef- 
ficients C„,: C = (1/M) J2m "i''^ ^ being the number of 
network mitochondria. 

FUNCTIONAL CONNECTIVITY 

To determine the "functional" connectedness of the mitochon- 
drial network, it is essential to determine the correlation coef- 
ficients between mitochondrial TMRE signals for each pair of 
mitochondria within the cardiomyocytes. For that matter, we 
followed a procedure by Eguiluz et al. (2005) to extract func- 
tional complex biological networks. This procedure abides by the 
recently introduced notion of the "spanning cluster" of oscillating 
mitochondria which stretches over the whole myocyte and does 
not necessarily involve all mitochondria (Aon et al, 2004, 2008a) 
as well as the concept of the dynamically changing major cluster 
of oscillating mitochondria described above (Kurz et al, 2010a). 
Two mitochondria were defined as being functionally connected 
at time t when the correlation coefficient of their TMRE signals in 
the time window around t was higher than a fixed cutoff value 
(see also Figure 3A in the results section). Consequently, cluster- 
ing coefficients were determined for every individual mitochon- 
drion and, for each considered myocyte (AT = 9), mean clustering 
coefficients C(t) (i.e., averaged over all mitochondria at time t as 
in Figure 3 or averaged over all major cluster mitochondria as in 
Figure 4) were determined for each point in time. 

RANDOM NETWORK 

For each myocyte, a random network was constructed that 
included the same number of mitochondria with the corre- 
sponding myocyte, using the Erdos-Renyi model (Erdos and 
Renyi, 1960). If Dm(f) presents the number of links of mito- 
chondrion m to its topological neighbors at time f (thus 
representing the dynamic degree of m), then each pair of 
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FIGURE 2 1 Time-dependent mitochondrial "mean field" 
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maps for 9 cardiac myocytes. The left hand side shows the correiation 


maps at the time point with the highest networl< mean correlation c(t) (see 


main text), whereas the right hand side shows those maps at the time 


point with the lowest mean correlation. Mitochondrial mean correlation 


values range between —0.4 and 0.8, the time-averaged mean correlation 


coefficient was found to be 0.43 ± 0.07. Also, some myocytes [myocytes 


(d) and (g)] show a prominent difference in overall correlation, while 


myocyte (i) does not show too pronounced correlation variability. 



mitochondria is connected with a dynamic probability p(t) = 
'^^D„,(t)/ {M(M — 1)). To allow comparison with the real 
mitochondrial network, the total number of connections of the 
random network at each time point remains the same as that of 
the real network. In such a network, most mitochondria have 
approximately the same number of links close to the average 
number of links at time f, ^,„D,„(f)/M (Barabasi and Oltvai, 
2004). 

STATISTICS 

Wavelet analysis, correlation analysis, and fitting routines were 
performed using Matlab v7. 6. 0.324 (R2008a). Further statistics 
were performed using OriginPro 8 SRO v8.0724 (B724). 

RESULTS 

INTER-MITOCHONDRIAL "MEAN FIELD " CORRELATION 

Mean-field coupling in complex non-linear biological networks 
of weakly coupled oscillators, such as a network of mito- 
chondrial oscillators (Aon et al, 2006) or in other systems of 
coupled chemical or physical oscillators, assumes a global, aU- 
to-all-coupling with a respective mean coupling constant (Kiss 
et al., 2002; Rougemont and Naef, 2006). In analogy to such a 
"mean field" approach, inter-mitochondrial correlation proper- 
ties were examined at time t through signal correlations c'™' (f ) 
of individual mitochondrion m with each of its topological 
network neighbors m; (Figure 2). The average correlation coeffi- 
cient ^iC^'^\t)/Nm for each mitochondrion m was then again 
averaged over all network constituents to determine the over- 
all dynamic inter-mitochondrial mean correlation c(t) as c{t) = 

For each myocyte, correlation maps were created to visualize 
the distribution of average correlation coefficients c-™'(f) within 

the myocyte at time t. To achieve this, values of c['"'(f) were 
placed at the pixel-positions of the respective mitochondrion and 
missing inter-mitochondrial pixels were interpolated using the 
"griddata" function in Matlab v7.6.0.324 (R2008a). In Figure 2, 
only the time points with maximum and minimum values of c(f ) 
for each respective myocyte are represented. It can be seen that, 
for some myocytes with maximum averaged correlation over all 
mitochondria, individual mitochondrial correlation values close 
to 80% can encompass almost the entire myocyte whereas the 
same can be valid for values of the correlation coefficient close 
to 0 (e.g., mitochondrion (d) and (g) in Figure 2). Averaged 
over all myocytes {N = 9) and time points, the mean correlation 
coefficient was found to be 0.43 ± 0.07. 

MITOCHONDRIAL NETWORK CLUSTERING 

The mean functional inter-mitochondrial clustering coefficient 
was determined for different correlation cutoffs (see Figure 3A). 
Naturally, high correlation coefficient cutoffs only involve few 
functionally connected mitochondria whereas lower cutoffs 
involve a higher number of mitochondria that are thought of 
as being connected. Therefore, the mean clustering coefficient 
increases for lower cutoffs. The ensuing fitted curve / (red 
line in Figure 3A, based on equation: /(f) = a ■ Exp [—t/b] + 
c) approximately crosses the 50% value of the mean clustering 



coefficients at correlation coefficients of 90% which also approx- 
imately corresponds to the median derivative value of the fitted 
curve. The fitting curve function parameters were determined as 
a = -7.81 ■ lO"*^ ± 1.33 ■ 10"*, b = -8.64 ■ 10"^ ± 0.13 • 10"^ 
and c = 77.65 ■ 10"^ ± 0.35 ■ 10"^. Accordingly, and throughout 
the paper, two mitochondria are thought of as being connected, if 
their correlation exceeds 90%. 

Time-normalization enables an estimate of the time evolution 
of the mean clustering coefficient (averaged over all mitochondria 
in the respective myocyte) for myocytes with unequal recording 
durations, and its comparison with that of random networks con- 
sisting of the same number of topological vertices and undirected 
links as the mitochondrial network of the respective myocyte 
(Figure 3B). Here, time-averaged mean clustering coefficients 
were found to be 0.500 ± 0.051 vs. 0.061 ± 0.020 for random 
networks, respectively. 

MAJOR CLUSTER PROPERTIES IN RELATION TO MEAN 
MITOCHONDRIAL CLUSTERING 

Individual mitochondrial signals are non-stationary in time 
and, therefore, wavelet transforms provide adequate means 
to examine the mitochondrion's signal temporal evolution 
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FIGURE 3 I Mean function clustering coefficient. (A) Time- and 
myocyte-averaged mean clustering coefficient for different cutoffs of the 
correlation coefficient as the defining threshold of functional 
connectedness. The red curve is an exponential fit curve of the form 
fit) = a ■ Exp [— t/t)] -I- c (for details and fit parameters see main text). (B) 
Time evolution of the mean clustering coefficient averaged over different 
myocytes (W = 9) (red line) and for random networks on the same set of 
vertices and maximum degree (blue line). Evidently, the mitochondrial 
network, as opposed to the randomly constructed network, shows 
significant functional clustering. 



(Grossmann et al., 1985). Following a recently described method- 
ology (Kurz et al., 2010a), mitochondria were sorted according 
to the dynamic behavior of their frequencies such that clusters 
of mitochondria with similar frequencies could be identified. 
Mitochondria from the major frequency cluster were subse- 
quently sampled and their mean temporal cluster area, amplitude 
and coherence determined (Kurz et al., 2010a). These spatio- 
temporal mitochondrial network properties are related to the 
functional topological network properties through the mean 
time-dependent clustering coefficient (see Figure 4). 

In Figure 4A, it is shown that the mean clustering coefficient 
increases with increasing mitochondrial cluster coherence. The 
latter takes values between 0 and 1 at each frequency and measures 
the degree to which two signals are synchronous or correlated: 
a coherence of 1 is equivalent to perfect synchronization. It is 
assumed that mitochondrial clusters possess a high dynamic sta- 
bility since their temporal coherence does not change significantly 
during the recording (Kurz et al., 2010a). The mean slope of the 
linear fit (red line in Figure 4A) was found to be 1.28 ± 0.09 
with the interception point at — 0.34 ± 0.06. This is not surprising 
since functional connectedness (as defined in the methods sec- 
tion) was chosen such that mitochondria are highly correlated. 
Higher degrees of clustering therefore correspond to a higher 



number of an individual mitochondrion's topological neighbors 
whose signals are highly correlated. 

Furthermore, mitochondrial cluster frequency and normal- 
ized cluster area (see methods section) were compared against 
mean clustering coefficients (Figures 4B,D): larger cluster areas 
show a higher degree of functional clustering due to the fact 
that an increase in mitochondria with similar frequencies natu- 
rally increases the number of mitochondria that are functionally 
connected. The mean increase of the mean clustering coefficient 
versus normalized mitochondrial cluster area (in the form of 
the linear fit in Figure 4D, red line) was found as 0.30 ± 0.01, 
the interception point at 0.28 ± 0.01. Also, the mean clustering 
coefficient was considered against the mean mitochondrial clus- 
ter frequency (Figure 4B) and it can be seen that clusters with 
increasing frequencies show lower clustering. The mean 

decrease of clustering versus major cluster frequency is —10.57 ■ 
10^^ ± 7.32 ■ 10^^%/mHz (the interception point of the linear 
fit at 0 mHz being 0.71 ± 0.02). The result is in agreement with 
prior results (Kurz et al., 2010a) showing that higher cluster 
frequencies coincide with smaller normalized cluster area, thus, 
with Figure 4D, lower mean mitochondrial clustering coefficients 
(compare Figure 4D with Figure 2 in Kurz et al., 2010a). 

Figure 4C presents the relation of topological clustering and 
the major cluster's amplitude, normalized to the maximum 
amplitude. Here, topological clustering reaches a plateau close 
to the value of the time-averaged mean clustering coefficient of 
50% for oscillation amplitudes that are lower than 80% of the 
maximum amplitude of the major cluster. But, for amplitude 
values close to the maximum value, clustering decreases to 
minimal values. The linear fit in Figure 4C has been achieved 
with/(f) = a ■ Exp [—t/b] + c (as in Figure 3A) — fit parameters 
were found as a = -52.01 ■ 10"** ± 26.97 ■ 10-^ fo = -7.29 ■ 
10^2 ± 0.28 ■ 10"^ and c = 51.10 • 10"^ ± 0.31 • 10"^ This is 
in agreement with prior observations (Kurz et al., 2010a) where 
it was shown that an increase in mitochondrial cluster amplitude 
corresponds to lower normalized cluster area, therefore, lower 
mean clustering coefficients (Figure 4 in Kurz et al., 2010a). 

DISCUSSION 

Mitochondrial networks in cardiac myocytes show evidence of 
interdependence between network structure and mitochondrial 
metabolism (Zorzano et al, 2010), as well as properties of a 
complex non-linear system under pathophysiological conditions 
(Aon, 2010). However, the only existing quantitative analysis of 
the mitochondrial network's topology is based on mitochondrial 
reticulum models (Sukhorukov et al., 2012), whereas other mod- 
els mostly considered the networks biochemical properties (Wu 
et al, 2007; Zhou et al, 2010). Thus, there has been no prior study 
of the networks functional connectedness through quantitative 
individual mitochondrial signal analysis. 

In the present work, we have shown that the mean correlation 
and the mean functional clustering coefficients in mitochon- 
drial networks are spatially and temporally variable. The latter 
were shown to be significantly higher than those for random 
mitochondrial networks with the same respective number of 
mitochondria and inter-mitochondrial topological connections. 
Furthermore, topological network properties were related to the 
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FIGURE 4 I Major Cluster Properties in Relation to IVIean Mitocliondrial 
Clustering. (A) Mean clustering coefficient vs. major cluster coherence. 
Though cluster coherence variability is limited, mean clustering still 
increases naturally for increasing cluster coherence. (B,D) Mean clustering 
coefficient vs. major cluster frequency and normalized cluster area. 
Clustering is shown to increase with cluster area and, therefore, number of 
cluster mitochondria. In analogy to Figure 2 in Kurz et al. (2010a), clustering 
then decreases with increasing cluster frequencies. (C) Mean clustering 



coefficient vs. normalized cluster amplitude. Clustering reaches a plateau for 
normalized amplitudes below about 4/5 of the maximum amplitude, 
whereas clustering decreases toward zero for maximum amplitudes. The 
latter effect might be explained by a more prolonged and incomplete 
repolarization in between oscillations when the relative size of the cluster 
increases whereas the plateau might signify a maximum of mitochondrial 
clustering. However, further amplitude decrease might be due to a failure of 
the mitochondrial network to reenergize. 



spatio-temporal properties of the major mitochondrial cluster 
with similar cluster coherence, frequency, normalized amplitude 
and cluster area. It could be shown that clustering increases with 
the percentage area of the cluster as well as the cluster's coher- 
ence. At first sight, this seems intuitively correct, however, the 
concept of clustering should not be confused with the notion 
of mitochondrial clusters in frequency. Nevertheless, increas- 
ing the number of mitochondria with similar frequencies, as is 
the case in synchronizing mitochondrial networks (Aon et al, 
2006), increases the number of mitochondria in the network 
that are functionally linked, and, therefore, mean mitochondrial 
clustering increases. Also, in agreement with previous results 
(Kurz et al., 2010a), mean mitochondrial clustering was shown 
to decrease with higher frequencies supporting the notion that 
smaller clusters of mitochondria oscillate with a higher frequency. 
It has been proposed that this process can be due to temporal 
limitations of diffusion-mediated inter-mitochondrial coupling 
that starts in small nuclei of excited mitochondria (Aon et al., 
2006) and progresses within the coupling medium (Kurz et al., 
2010a), therefore leading to a decrease in the mean cluster fre- 
quency through an increase in depolarization time, and/or due 
to coupled mitochondrial oscillators adjusting to a common 
oscillatory mode with larger clusters taking longer to synchronize. 



Additionally, it has been demonstrated that increased normal- 
ized cluster amplitudes are related to a decrease in mean mito- 
chondrial clustering. This less intuitive result can be explained 
by the fact that mitochondria in the major cluster show a more 
prolonged as well as incomplete repolarization in between oscil- 
lations when the relative size of the cluster increases. Then again, 
the increase of cluster-incorporated mitochondria leads to an 
increase of topological clustering (cf Aon et al., 2008b). Yet, 
Figure 4C indicates that this effect comes to a halt for ampli- 
tudes below approximately 4/5 of the maximum amplitude. In 
this amplitude realm, mean mitochondrial clustering remains 
constant at its time-averaged value for decreasing cluster ampli- 
tudes suggesting a stable maximum of topological connectedness. 
However, it should be noted that the effect of fluorescence inten- 
sity loss can be due to limitations of the TMRE readout that 
might be inadequate to track correlations in the state when 
redox equivalents supply is insufficient to reenergize the mem- 
brane. Under these conditions, energetic/redox impairment may 
eventually lead to cell death. 

Mitochondrial networks consist of regularly organized lattice- 
like mitochondria that are seen as metabolic network hubs with 
multiple connections to anabolic and catabolic cellular pathways 
(Aon et al., 2008b; Aon and Cortassa, 2012). Their coordination 
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was first dynamically examined in cardiac myocytes that were 
subjected to metabolic stress (Romashko et al., 1998; Aon et al., 
2003) where it was presumed that ROS-induced ROS release 
acts as a coupling mediator between mitochondria to coordi- 
nate local and global cell-wide mitochondrial network behavior 
(Zorov et al., 2000; Aon et al, 2004). The concept of scale-fi-ee 
networks and its application to the mitochondrial network's fi^nc- 
tional connectivity has been the focus of recent studies (Aon 
et al, 2004, 2006, 2007, 2008b; Barabasi and Oltvai, 2004), which 
support the idea that the mitochondrial network architecture is 
connected to its function (Viola et al., 2009) and that metabolic 
supply changes can be due to morphological changes of the 
mitochondrial network structure, and, specifically, the morphol- 
ogy of mitochondrial membranes (Dimmer and Scorrano, 2006; 
Mcbride et al, 2006). 

The present study also supports this concept in showing that 
topological clustering in mitochondrial networks is significant 
and relates directly to the network's spatio-temporal organiza- 
tion. It has also been shown that a mitochondrion's localization 
within the mitochondrial network has an effect on its func- 
tional properties (Lesnefsky et al., 1997; Kuznetsov et al, 1998). 
This is in agreement with the observation of strongly variable 
areas of mean-field correlation in mitochondrial networks as 
well as time-variant mitochondrial clustering. In addition, several 
studies could demonstrate that some individual mitochondria 
show independent responses to mitochondrial network dynam- 
ics (Loew et al, 1993; Duchen et al., 1998), again emphasizing 
mitochondrial functional heterogeneity and the complexity of 
the network's functional connectedness. The principle of func- 
tional connectedness applied in this study has also been used to 
describe brain functional networks that were shown to exhibit 
significantly larger clustering coefficients than those of random 
networks (Eguiluz et al., 2005); they are subject to intensive 
clinical research (see BuUmore and Sporns, 2009 for a review). 
Interestingly, correlations between spontaneous changes in brain 
activity generate very robust functional networks on timescales of 
seconds to minutes, i.e., in the same range as mitochondrial net- 
work dynamics (Greicius et al, 2003; Fox et al., 2005). Continuing 
this comparison, neuronal phase synchronization between cou- 
pled neurons was found to be highly dependent on clustering 
(Percha et al., 2005) and to possess self-organized critical dynamic 
properties (Siri et al., 2007). 

The present analysis is limited to a functional analysis of 
the mitochondrial network. Furthermore, the present study only 
accounts for spatio-temporal properties of mitochondria in the 
major frequency clusters but neglects smaller clusters. This can be 
justified since the major cluster usually comprises the majority 
of network mitochondria that are above a percolation thresh- 
old of approximately 60% for lattice-like organized networks (cf 
Aon et al., 2004). Nevertheless, the effect of long-range connec- 
tions between different clusters and, therefore, the modularity of 
the network, will be the subject of further studies. Methods that 
examine the effective connectivity of the system to evaluate the 
causal influence of each system element on the network's behav- 
ior (see Friston, 2011 for a review), could be used to assess the 
connectivity between mitochondrial clusters within the myocyte 
or between neighboring myocytes, and consequently to examine 



scaling effects of mitochondrial clustering. It should also be noted 
that only two-dimensional slices of the cardiac myocytes have 
been examined here, that consequently contained only parts of 
the three-dimensionally arranged mitochondria. However, we 
would expect lateral coordination among rows of mitochondria 
above and below those in the focal plane to be similar to the 
communication between adjacent rows in the 2D plane, since 
the mitochondrial network has a symmetrical quasi-square lat- 
tice organization. The main source of anisotropy is the lateral 
versus longitudinal orientation of the mitochondria arrangement 
between the myofilaments, as evidenced by preferential depolar- 
ization propagation in the mitochondrial network along the axis 
parallel to the cardiac myocyte myofilaments (Kurz et al, 2010a). 

In summary, the results of this study indicate that cardiac 
mitochondria constitute a collection of coupled network com- 
ponents that are subordinate to dynamic changes in their non- 
random functional connectedness. 
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